Nonexponential motional damping of impurity atoms in Bose-Einstein condensates 
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We demonstrate that the damping of the motion of an impurity atom injected at a supercritical 
velocity into a Bose-Einstein condensate can exhibit appreciable deviations from the exponential 
law on time scales of 10 -5 s. 
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The decay of an unstable quantum system has long 
been known to deviate from exponential law for both 
very short and very long times [1] - [4]. The short-time 
deviation from the exponential decay gives rise to cither 
slowdown or speedup of the decay by frequently repeated 
measurements, known, respectively, as the quantum Zeno 
effect (QZE) [1] and the anti-Zeno effect (AZE) [2,3]. A 
general unified theory of the QZE and AZE has been 
given in Ref. [4]. Among the proposed realizations of 
the QZE/ AZE are photodetachement of negative ions 
[5] and radiative decay of excited atoms in cavities [3], 
photonic band gap structures [6] or in the presense of a 
Bose-Einstein condensate (BEC) [7]. The first unambigu- 
ous observation of nonexponential decay in an unstable 
quantum system has been reported in [8] , followed by the 
demonstration of the QZE and AZE [9]. These experi- 
ments have measured the escape of cold atoms from wells 
of an accelerating periodic potential induced by a stand- 
ing light wave with varying frequency. Qualitatively, es- 
cape from a trapping potential resembles nuclear alpha- 
decay [10] rather than decay via quanta emission into a 
bath, as in radiative decay [11] or in beta-decay [12]. 

The difficulty impeding the demonstration of nonex- 
ponential decay via quanta emission has been the short 
non-Markovian (memory) time of the relevant (electro- 
magnetic or leptonic) bath, which is the time-scale of 
the effect [4]. Here we propose the emission of sound 
quanta (phonons) by atoms into a BEC as a candidate 
process for the observation of nonexponential relaxation, 
taking advantage of the rather long memory time of the 
condensate. The detailed understanding of such deco- 
herence processes is important for the envisaged use of 
atomic BECs in metrology and interferometry [13]. 

Let us consider an "impurity" atom moving in an 
atomic BEC. Atoms of another isotope or the same iso- 
tope but in a different internal (hyperfine) state can be 
viewed as impurities as long as their density is small 
enough not to modify considerably the BEC excitation 
spectrum. At "supercritical" velocities, namely, above 
the speed of sound in the BEC, the impurity atom is de- 
celerated due to phonon creation in the BEC. The rate 
of such a process according to the standard Fermi golden 
rule (i.e., assuming exponential decay of the amplitude of 
the initial state) has been calculated for both a uniform 



BEC [14] and a harmonically trapped BEC [15] and has 
been determined experimentally [16]. We shall step be- 
yond this approach and calculate more generally the time 
evolution of the impurity-atom motion in such a system. 

The state vector of the system can be written as (we 
set h = 1) 
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where the initial state |in) corresponds to an impurity 
atom of mass mi moving at the velocity Vi in a BEC 
with no elementary excitation, and |q) denotes the state 
where one elementary excitation of the BEC with the 
momentum q is present (correspondingly, the impurity 
momentum is changed to m{V \ — q). The initial condi- 
tions are, naturally, 



t (0) = l, 0,(0) =0. 
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The set of equations describing motional damping of the 
impurity atom is 
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Here n is the uniform BEC density, gi2 is the effective 
interspecies coupling constant, u q = {[E h f (q) / t(q) + 
lJ/2} 1 / 2 and v q — (u 2 — l) 1 / 2 are the Bogoliubov trans- 
formation coefficients, and e(q) — [Ejj F (q) — /J 2 ] 1 ^ 2 is 
the energy of the elementary excitation with momentum 
q [17]. Here we have introduced the Hartree-Fock exci- 
tation energy Ehf = ? 2 /(2 m 2) + /i and the chemical 
potential of the BEC /i = Anazzn/mi, a-ii being the 
intraspecies s-wave scattering length for the condensed 
atoms of mass to 2 • Correspondingly, the speed of sound 
in the BEC is c s = •\// i / m 2- 
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To proceed, we have to reckon with the renormalization 
of the interspecies coupling constant [18]. The renormal- 
ized constant g\2 is expressed in terms of the bare cou- 
pling constant 312 as 



312 = 9i2 



1 + 2m 5l2 (27r)' 3 / d 3 qq 
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with gi2 — 27rai2/m, ai2 being the interspecies s-wave 
scattering length and m = miiJi2/(mi + 7712) being the 
reduced mass. In the approximate (perturbative) solu- 
tion of Eqs. (3, 4) we shall keep the terms up to the 
second order in the bare constant 312- 

The easiest way to solve Eqs. (3, 4) is by the Laplace 
transformation (cf. Ref. [6]). We adopt the interaction 
representation, wherein the probability amplitude of the 
initial state is 



a(t) = a i „(i)exp[i(miV 1 2 /2 + .gi2n)i]. 
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The algebraic solution for the Laplace transform of a(t), 
a(s) = J °° dt exp(— st)a(t) has the form 



y(s) = [s + n(«)] 1 
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A(q) = e(q)+q 2 /(2mi)~qV 1 being the energy mismatch 
between the the states |in) and |q). The second term in 
the square brackets in Eq. (8) arises from the coupling- 
constant renormalization in Eq.(5) and compensates for 
the ultraviolet divergence of the first term. This com- 
pensation is completely analogous to that of the electron 
mass renormalization in calculations of the radiative shift 
of an atomic optical transition [19]. 

Equations (3, 4) yield over a broad time interval, ex- 
cluding very short times, exponential decay of a(t) oc 
exp[— (7/2 + iw s )t], with the rate 7 and the frequency 
shift uj s 



7 = lim 2Rcft(s) 



lu s = lim Im O(s). 

s— »0 
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The relaxation rate 7 can be calculated within the exact 
Bogoliubov theory using Fermi's golden rule [14,16]. 

It is possible to obtain a nonexponential analytical so- 
lution for a(t) in the Hartree-Fock (HF) approximation 
[20] (hereafter we label all the quantities in this approxi- 
mation by HF): e H F(q) = E HF (q), u qHF = 1, v qHF = 0. 
Then the integral in Eq.(8) can be evaluated analyti- 
cally, yielding in Eq. (9) the exponential decay rate 



1HF 



7o [1 



"HF 



if Vl > c HF = 7277/^ and 



zero otherwise. Here 70 = Aira^nVi is the collison rate 
calculated for the impurity atom using the Fermi golden 
rule in the limit 022 — > (an ideal BEC). Finally, in the 
case V\> Chf we obtain 
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where K = (mV 2 /2 — h)/^hf is the energy of collision 
of the impurity atom with a condensate atom in their 
center-of-mass frame (including the correction due to the 
BEC chemical potential jit) scaled to jhf- By inverting 
such a Laplace transform [21] we get 
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£i¥>(£iVt) - 6^(6^)1 , (11) 



where <p(z) — exp(2 2 )erfc(— z), erfc(z) being the comple- 



mentary error function, and £i )2 = — S/2 ± ^/S 2 /4 — iA 
are the roots of the quadratic equation £ 2 + S£ + iA = 
0. The coeffi cients of the latter equation are S = 
7o V / V( 2to ^i 2 ) and A = mV?/2 - n. 
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FIG. 1. Numerically calculated logarithm of survival prob- 
ability P(t) of the initial state of impurity atoms in a S7 Rb 
BEC plotted versus dimensionless time -yt for n — 10 14 cm -3 
(c s = 0.2 cm/s); mi = rri2, (I12 = 3a22- Filled circles: 
Vi = 3c s , 7 = 1.1 • 10 3 s -1 . Open circles: Vi = 7c s , 
7 = 4.4- 10 3 s _1 . Solid line: exponential law exp(— ft). Inset: 
typical behavior of logarithm of P(t) for the case of subcriti- 
cal Vi (practically any Vi < 0.7c s ). Note the different scaling 
(by fi) of the horizontal axis of the inset plot. 

Although our numerical results (Fig. 1) show that the 
HF approximated solution Eq. (11) is rather crude, it 
nonetheless provides a qualitative guidance to the phys- 
ical behavior. Simple scaling considerations lead us to 
the conclusion that oihf depends on two parameters: 
the dimensionless time variable "fHFt and K [Eq. (10)]. 
Equation (11) predicts that at t — > the decay is 
more rapid than exponential, so the survival probabil- 
ity P(t) = \a(t)\ 2 behaves in the HF approximation as 
PnF(t) ~ 1 — 4Re(£i + 6)\A/ 7r i becoming exponen- 
tial at larger times, Ph f (i ) ~ exp(— fHFt). The devia- 
tion from exponential decay is appreciable only for K < 1. 
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There are two ways to attain K < 1. One is to take a 
small difference between the impurity atom velocity V± 
and the critical velocity, but this would reduce the damp- 
ing rate, which may be experimentally inconvenient. A 
much better way is to strive for a large interspecies scat- 
tering length ai2, as discussed below. 

The results of our numerical calculations based on 
Eqs. (7, 8), which use the exact expressions for e(q), u q 
and v q instead of the HF approximation, clearly reveal a 
deviation from exponential decay for small times. Under 
such conditions, frequent measurements would accelerate 
the decay, causing the anti-Zeno effect (AZE). Alterna- 
tively, one may accelerate the decay by periodically mod- 
ulating the coupling of the initial state to the continuum 
[4], instead of repeated projective measurements. This 
can be done by changing the impurity velocity using a 
sequence of Bragg or Raman laser pulses [22]. 

The inset to Fig. 1 shows that if V\ < c s , the sur- 
vival probability first decreases and then approaches the 
constant value of about 0.85, almost independent on V\. 
This behavior reveals the physical reason for the short- 
time nonexponential decay: the initial conditions Eq.(2) 
imply that, initially, the impurity atom is surrounded by 
no virtual phonons, while in the steady state, the im- 
purity atom must be surrounded by a cloud of virtual 
phonons (cf. the polaronic effect for electrons in a crys- 
tal [23]). Thus the nonexponential stage of the decay is 
associated with the formation of such a phonon cloud. 

The faster the impurity atom moves, the weaker its 
coupling to the phonon cloud is. Therefore the decrease 
of lu s corresponds to the vanishing of nonexponential de- 
cay effects as V\ increases. This behavior is displayed in 
Fig. 2 by the numerically calculated [from Eq. (9)] decay 
rate 7 and frequency shift u s . 

Our numerical studies of Eqs. (3, 4) always yield de- 
cay acceleration at short times, typically on the scale of 
ICE 5 s. But should one not expect, from general consid- 
erations [1], P(t) = 1 — const • t 2 at t — > 0, in accordance 
with the QZE? To answer this question, we should ap- 
ply the general theory developed in Ref. [4], whereby 
the short-time behavior is determined by the spectrum 
(i.e., the dependence on the emitted quantum energy e) 
of the reservoir response G(e). This spectrum is given by 
the interaction matrix element squared multiplied by the 
density of the reservoir states. In our case, we find that 
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monotonously increases with the emitted phonon energy 
e. According to Ref. [4], if the energy uncertainty ~ t^ 1 
associated with the finite observation time t covers the 
energy range where G(e) increases, then decay accelera- 
tion (AZE) takes place. However, our approach [Eqs. (3, 
4)] leading to Eq.(12) is valid only for small transferred 
momenta. If q^r$ , where ro ~ ICE 7 cm is the char- 
acteristic radius of the interatomic potential, we cannot 
consider a\2 as constant any more. Instead, the interac- 



tion matrix element decreases with q in this range. In the 
inset of Fig. 2 we schematically display the spectrum of 
G(e), including its decreasing part, whose detailed calcu- 
lation is beyond the scope of this Letter. This spectrum 
implies that at very short times (< 10~ 9 s), the energy 
uncertainty broadening covers the whole profile of G(e), 
thereby giving rise to the QZE. 




FIG. 2. Numerically calculated frequency shift u) s and the 
exponential decay rate 7 (scaled to 70) versus the impu- 
rity dimensionless velocity Vi/c s . Both LUs/'Jo and 7/70 dis- 
play a universal behaviour, independent of the BEC density 
and the atomic species. Inset: schematic representation of 
the spectrum G(e); solid line: the increasing part given by 
Eq.(12); dashed line: the remaining part, which decreases at 
e>l/{mrl). The arrows indicate ranges of the energy uncer- 
tainty corresponding to AZE and QZE. 

The parameters used in Fig. 1 for a BEC of 87 Rb may 
correspond to impurity atoms of the same isotope but 
in a different hyperfine state, obtained by a short Raman 
pulse. This method of impurity atom admixing conforms 
to the initial conditions of Eq.(2). However, to reach ap- 
preciable nonexponentiality by this method, one has to 
enhance interspecies scattering either by means of inter- 
species Feshbach resonance or via laser-induced dipolc- 
dipole interactions [24]. Another possibility is the use of 
two-isotope mixture, for example, a BEC of 87 Rb atoms 
with admixed fermionic 40 K atoms. Such a choice is of 
particular interest, since the large interspecies scattering 
length in this mixture [25], as well as the mass ratio be- 
tween these two elements seem very promising for exper- 
imental search of nonexponential decay effects. We note 
that the condition on the impurity velocity in this case is 
opposite to that of the aforementioned case of impurities 
generated from a BEC by a Raman pulse. Indeed, the 
40 K atoms at rest co-exist with the 87 Rb BEC for a time 
long enough to form virtual phonon clouds around them. 
If then the 40 K atoms suddenly acquire, by the action 
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of a Bragg or Raman pulse, the velocity V\ 3> c s , the 
initial conditions, instead of Eq.(2), should assume the 
pre-existence of the phonon cloud (prior to the pulse), 
which vanishes when the impurities attain high (super- 
critical) velocities. Thus we expect that in the case of 
a two-element ultracold mixture, noncxponential decay 
features are most pronounced for V\ » c s . 

To conclude, we have outlined the possibility of ob- 
serving deviations from exponential decay for unstable 
momentum states of impurity atoms moving in a BEC 
with a supercritical velocity on time scales of 10~ 5 s. The 
effects of finite temperature and bosonic enhancement in 
this process will be a subject of a separate work. 
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